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We find the shock wave solutions in a class of cosmological backgrounds with a null 
singularity, each of these backgrounds admits a matrix description. A shock wave solution 
breaks all supersymmetry meanwhile indicates that the interaction between two static D0- 
branes cancel, thus provides basic evidence for the matrix description. The probe action of 
a DO-brane in the background of another suggests that the usual perturbative expansion 
of matrix model breaks down. 
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It is an important open problem to formulate string/M theory in a time-dependent 
background in general, and to study what happens to cosmological singularity in particular. 
The recently proposed matrix description of a cosmological background with a null 
singularity appears to be an advance forward, this construction was generalized to a class 
of backgrounds in 0] . An earlier example of cosmology with a null singularity is proposed 
in HU and is later studied by many authors j5|-0], they find that this singularity is highly 
unstable due to gravitational back-reaction. It seems to us that models studied in J^,0 
and in the present note are devoid of such problems. 

A cosmological background with a null singularity is not necessarily stable, even when 
it admits a perturbative string description P-^]. It was found that a smeared shock wave 
changes the nature of the singularity, this happens in 2 + 1 dimensions. Our solution 
presented below is also smeared in one or more directions, but localized in other directions, 
there is no sign of instability. 

Matrix theory |8|,[| description predicts that there is no interaction between two static 
DO-branes, this is verified by our solution. The matrix action does not possesses linear 
supersymmetry (there is always non-linearly realized supersymmetry corresponding to shift 
of fermions), we show that the shock wave solution breaks all supersymmetry, in accordance 
with the matrix action prediction. Does the fact that there is no net force between two 
static DO-branes imply that we can expand the effective action in terms of relative velocity? 
The probe action in the background of a shock wave shows that this is not the case. This 
has far-reaching consequence: the usual perturbative expansion in matrix model breaks 
down. 

The background studied in [Q is a flat 10D space-time with a linear dilaton propor- 
tional to the light-cone coordinate x + . Its 11 dimensional M theory matric reads 

ds 2 u = e 2 ^ + / 3 ds 2 10 + e-^ x+ /\dY) 2 . (1) 

The 10 dimensional part shrinks when x + — > — oo, thus there is a null singularity a finite 
distance away in terms of the affine parameter. It is argued in |l] that the matrix model 
works for this background with a time dependent action, where the matrix time t is chosen 
to coincide with x + . 

The above solution can be generalized to the following class 

ds 2 = e 2ax+ (-2dx+dx- + (dx l ) 2 )+e 2(3x+ (dx a ) 2 

(2) 

i = l,...,9-d, a = 9-d+l, ...,9, 
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where the transverse directions are grouped into two sets. There are other generalizations 
in 0, but here we restrict our discussion to the above two cases. 

It is believed that for each of these backgrounds there is a matrix description, which 
may provide a theory for resolving the singularity. For this claim to be valid, it is important 
to check whether these backgrounds are stable. Here we take an approach similar to that 
of ||. Consider there be some gravitational sources, localized in x~ dimension as well as 
in the transverse dimensions. The existence of such sources will change the background. 
We want to see how large the change can be. First we consider the easiest case, with 
the background metric given by (|IJ): It is convenient to treat the x~ and Y dimension as 
periodically identified, i.e. 

x~ ~ x~ + 2nR, Y ~ Y + 2nR' 
Consider the perturbation represented by stress tensor 

T__ = Y,P-a e~ 2Qx+ 6(x ± - x ±a ), (3) 

a 

where the x + dependence is introduced by energy conservation (this is for the background 
(P; for (0), the x + dependent factor is exp[— ((9 — d)a + d(3)x + }), and x± represents all 
the uncompactified transverse dimensions. The above stress tensor is due to multi-sources, 
located at x± a in the transverse directions. The stress tensor is smeared in the x~ and 
Y directions, thus p- a is not really the longitudinal momentum of the a-th particle, the 
longitudinal momentum should be p- a 27rR2nR' , where R is the radius of x~ and R' the 
radius of Y. We make the ansatz that the metric takes the following form 

ds 2 = -2e 2ax+ dx+dx~ - 2e ex+ f(x~ ,x A )(dx-) 2 + e 2aAX+ (dx A ) 2 , (4) 

where A = i, a in @, 6 is a constant to be determined. This ansatz can solve Einstein 
equations with source (|3D, and requires more non- vanishing components in stress tensor if 
the source is not smeared in direction Y (or in directions x a ), as we shall see shortly. The 

component reflects the effect the source has on the geometry. If g vanishes, we 

get back to the original metric we started with, which should satisfy the vacuum Einstein 
equations with the constraint on the parameters: 

cla{oia — 2a) = 0. (5) 

A 
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For metric (|l]), we have 



a = ai = Q/3, i = l... 8; a 9 = -2Q/3. (6) 
In the case of (||), there are 

ai = ct, a a = (1 ± 3/Vd)a. (7) 

To find a solution with a source (|3|), we use the orthonormal basis 

e~ = dx~ 

e + = e 2ax+ dx+ + e 6x+ dx~ (8) 
e A = e *AX + dx i_ 

Up to symmetry, the non-vanishing components of spin connection are 

U++ = 9e i9 - 2a ^ x+ dx~ 

u+a = e {B - aA)x+ d A fdx- - a A e (aA+e - 2a)x+ fdx + (9) 
u A + = a A e {aA ~ 2a)x+ dx A . 

The non- vanishing curvature 2-forms are 

0++ = 9(6 - 2a)e {e ~ Aa)x+ fe + A e~ + (6 — a A )e {e - 2a ~ aA)x+ d A fe l A e~ , 

Q+ A = [e^ e - aA -^ x+ d A d B f - a A (a A - 2a)e ( - 29 -^ x+ f 2 5 AB + a A e^ e - 2 ^ x+ d_f 8 AB ) 

e B A e~ + (9 - a A )e (e - 2a - aA)x+ d A fe + A e~ + a A e {e - 2a ~ aB)x+ d B fe A A e B 

+ a A (a A + 6- 2a)e (e - 4a)x+ fe A A e+ 

Q A + = -a A (a A - 2a) e~ Aax+ f e+ A e~ + a A (6 + a A - 2a)e^~ 4a)x+ f e A A e~, 

n A B = a A e ^ aA - 2a > + e A Ae~ — a B e^ aB ~ 2 ^ x \ B Ae~ — 2a A a B e^~^ x+ e A A e B . 

(10) 
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In the Lorentz frame, the non-vanishing components of the Ricci tensor are 
R++ = -Y,aA(a A -2a)e- 4ax \ 

i 

R— = - a A {a A - 2a)e i - 20 -^ x+ f 

A 

+ ]T e^~ 2a ^ x+ d A d A f + a A e^ e - 2aA > + cL/, 

A A 

R+- = [6(6-2a) + J2 a A(0 + a A -2a)] e ^~^ x+ (11) 

A 

Ra+ = o, 

R A _ = (6-2a A + J2 «b) e( e ~ 2a -^ x+ d A f, 

B 

Rab = ~2a A (9 + Y, a c~ 2a)e ( *" 4a)x+ / S AB . 
c 

Because the transformation from the orthonarmal basis to the coordinate basis does not 
mix the spatial spatial components, we see from flTT| ) that the condition T A b = implies 

6 -2a + ^a A = 0. (12) 

A 

This condition also implies that R A - = if a A = a. This is the case for the background 
(0) when A = i. If a A ^ a, we will have to postulate d A f = 0, namely we need to smear 
the source in these directions. If we do not, a non-vanishing component T A _ is required, 
this is just fine, since R A - thus T A _ is a total derivative, the component of momentum 

Pa ~ / T A -yfgZdx~d 9 x (13) 



automatically vanishes. However, we will assume that as long as a A ^ a, the source is 
smeared in direction x A . Finally, Taking (|12|) as well as (|5|) into account, and changing 
back to the coordinate basis, we find the only non- vanishing component of Riccci tensor 
R__ 

R— = e {B ~ 2aA)x+ d A d A f + a A e {B ~ 2aA)x+ d-f. (14) 

A A 

From now on we focus on the two cases ([!]) and (||) only. In the case of ([!]), / is 
independent of a; - and Y, and 

9 = -4Q/3. (15) 



The remaining Einstein equation is 

8e - 2 Qx+ A/ = Ki 2 J2p- a e- 2 ^ + 5(x ± - x ±a ), (16) 

a 

where A is the Laplace operator in the 8 transverse dimensions. Finally we get the solution 

/ = -E l5 r lP " a - i6 > ( 17 ) 

^ 48S7|Xj_ — X± a \ 

where sr denotes the volume of the sphere S 7 . 

The surface defined by x + = constant becomes spacelike now. So we will get a spacelike 

singularity. Quantitatively, when g__ and become the same order, the geometry is 

severely changed. Consider only one source, there is a radius: 

r = e -°- + (i^L)i (18) 

KllP- 

The condition r ^> tq defines the asymptotical region. When x + — > oo,ro — > oo. It 
seems dangerous. But there is some subtlety here. First, When x + is finite, r$ is also 
finite, and we do not need to worry about the change as long as we are far away from 
the source. As x + grows, r$ vanishes, and so the future is changed little. Second, if we 
want to detect the change of background, we must introduce another observer. Here we 
put another graviton into the spacetime. The change of background is reflected by the 
change of action of the second graviton. As we will see later, this picture breaks down as 
x + — » — oo. So even the nature of the singularity is changed, we cannot detect the change 
on the classical level. To further understand it, perhaps we should go to matrix theory. It 
is important to ask whether matrix theory is right here. In matrix theory, two gravitons 
correspond to block diagonal parts of a matrix. They interact via the off diagonal part. 
Spacetime emerges as we integrate out the off diagonal part. When x + — > — oo, we reach 
the singularity in sugra, and a week coupling super Yang-Mills theory in matrix. As argued 
by [0, simply integrating out is not allowed. We must deal with the whole non-Abelian 
degrees of freedom. To us, there is no sign of breakdown of matrix theory. More rigorous 
matrix calculation is needed to say anything more. 

In a flat background as well as in a anti-de Sitter background, stability of a perturba- 
tion amounts to requiring that the space-time remains asymptotically flat or anti-de Sitter, 
since definition of S-matrix or boundary correlators demands this. In a time-dependent 
background with a big-bang singularity, intuitively, the space-time structure is required to 
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remain the same where there is no singularity as time approaches infinity, since we may 
invoke perturbative string picture there. However, at the point where there is singularity, 
the definition of observables is not clear yet, so we do not know how to define stability. 
Nevertheless, the kind of divergence found in || is not allowed, happily that does not occur 
here. 

Unlike the case discussed in 0, the change of the geometry is localized here. In the 
region far away from the sources, space-time is almost the same as the original background. 
In this sense, the background is stable. In the case of (0), the calculation is similar. 
Provided that we compactify x~ and x a and smear the source in these directions, all 

components of Ricci tensor vanish except R The solution of / is also similar, r 6 in 

the denominator in (|17|) is to be replaced by r 7 ~ d . If we take the radius of x~ to infinity, 
the stress tensor should also be localized in this direction, the source will contain a factor 

5(x~), so does the component g of the metric. S(x~) is also localized, it will not trigger 

instability. 

What we have obtained is noting other than the smeared shock wave representing 
the classical background of multi DO-branes. Eq.(^7|) clearly indicates that there is no 
net interaction among DO-branes. We now discuss whether this background still preserves 
supersymmetry. 
• Supersymmetry 

In 11 dimensions, the only SUSY transformation of interest is that for the gravitino 
Si/^n = D^e. The components of spin connection of interest are o>_ and 

uo- = 29e {9 - 2a)x+ fT+- + e {e ~ aA)x+ d A f T A ~, 

(19) 

u A = 2a A e (aA ~ 2a)x+ T A+ - 2a A e^ +aA ~ 2a)x+ f T A ~ . 
Then the condition D+e = says that e is independent of x + . For the condition 

D_e = (d- + -o;_)e = (20) 

multiplied by T~, it becomes 

(8- + -ee i9 - 2 ^ x+ f)T~e = (21) 

The second term is dependent of x + , while the first term not. Thus they both vanish, 
implying that e is also independent of x~ , and that e is subjected to the constraint T_e = 0. 
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Similarly, multiply -D^e = by T_, we end up with T^e = for all transverse directions, e 
has to vanish. Hence we have shown that the existence of sources breaks all supersymmetry. 
• Interaction 

We now switch to the convention in which x~ = x — t, hence the momentum p_ 
is always positive, and / is always negative. Consider the case that there is only one 
source. The solution looks much the same as the Aichelburg-Sexl metric. Hence we follow 
the approach of |TD[ to discuss the probe action of graviton of momentum N%/R in the 
background of graviton of momentum N\/R. When N\ is large, the first graviton is treated 
as a classical source and produces the metric we have got. Then take the other particle as 
a probe particle. To get the action of the probe graviton in this background, the trick is 
to begin with the action for a massive scalar in eleven dimensions 



S = -m J dri-G^afx")* 

= -rnj dr(-2e 2ax+ x- -^e 2aAX+ (x A ) 2 + 2e dx+ fx~ x') 1 * . 



(22) 



A 

A dot denotes d T . Carry out a Legende transformation on x~: 

p_ = m(e 2ax+ - 2e Bx+ fx~)(2e 2ax+ x- - e 2aAX+ v 2 A + 2e dx+ fx~ x~)~^ . (23) 
In our case p- is to be fixed. Taking the limit m — > implies that Gu iV x il "x 1 ' — > 0. Then 



= " 777~0x~+ " ( 24 ) 



e 2ax - \je^ + +2J2 A e( 9 + 2a ^ x+ (x A ) 2 
The appropriate Lagrangian is the Routhian 



£'(p_) = -n(p-) = C-p-x~{p-). (25) 

When m — > with fixed p_, 
C'(p-) — > —p-X~ 



p_{-Y^ e 2 ^~^ x+ (x A ) 2 - - e - &a+e [J2 e 2aAX+ (x A ) 2 } 2 f ( -> 6 ) 

2 A A 

+ 0(f 2 x\)}. 



The interaction terms all involve relative velocity v. If the source and the probe stay 
relatively static in the transverse space, there is no net interaction. This has to be the case 
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if matrix model is correct, since two static DO-branes form a solution in matrix model, and 
induces no interaction terms in the matrix model action. 
When the parameters take the values in (||), we get 

£'(p_) - P-{\^^f + e~ 2Qx V) 2 ] 

i 

- zE e ~ Qx+ ^ l ) 2 + e ~ 3Qx+ ( i9 ) 2 ] 2 f (27) 

i 

+ 0(fx\)}. 

We see that as x + — > — oo, those "higher order" terms become larger and larger, it is a 
bad idea to do expansion in terms of the transverse velocity. 

The double expansion in v and r _1 corresponds to perturbative expansion in the 
matrix model, the breaking-down of this double expansion implies that even when one 
is treating interaction between two DO-branes, one should no do perturbative calculation 
in matrix model. It is rather curious that this breaking-down happens, since the matrix 
model or the matrix string theory in question becomes weakly coupled at early times. 
Our interpretation of this breaking-down is the following. Although the theory becomes 
weakly-coupled, the picture of two separated DO-branes becomes worse and worse at early 
times, since the off-diagonal modes become lighter and lighter, this can be seen for instance 
from the interaction term 

Rg-hi[X\ X J } 2 = Re 2Qt tr[X\ X J ] 2 , (28) 

where t = x + . The mass of the off-diagonal modes is proportional to e^*. The breaking- 
down of velocity expansion is simply due to appearance of new light modes. Thus, although 
the DO-brane "moduli" (position) description breaks down, the matrix model does not, 
one simply has to deal with all non-abelian degrees of freedom. This is how matrix model 
resolves the singularity. 
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